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i  coupled  oceanic-atmospheric  boundary  lay er>  model  which 
provides  single- station  prediction  capability  *Cs  evaluated 
relative  to  wind  stress  and  oceanic  mixed  layer  depth  obser¬ 
vations.  The  model  is  initialized  and  verified  using  data 
obtained  during  the  1983  nixed  Layer  Dynamics  Experiment 
(MI1DEX).  Model  computation  of  friction  velocity  (u*)  and 
oceanic  mixed  layer  depth  are  compared  with  observations 
from  both  atmospheric  frontal  and  non-frontal  synoptic  situ¬ 
ations.  Favorable  results  of  model  u*  predictions  are 
achieved  although  in  some  cases  they  are  slightly  higher 
than  observed  turbulence  derived  values  of  u*.  Mixed  layer 
predictions  are  very  close  to  the  observed  except  in  frontal 
regions.  Proper  ABt/OBL  model  initialization  of  boundary 
layer  values  is  critical  in  order  to  achieve  favorable 
results.  K  ....  :  y 
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I.  INTRODUCTION 


Naval  oceanographers  and  meterologists  obtain  a  great 
variety  cf  information  to  help  them  in  making  atmospheric 
and  oceanic  forecasts.  Extensive  use  of  large  shore  based 
computers,  remote  sensing  equipment  and  world-wide  communi¬ 
cations  networks  make  the  forecaster  rely  on  remotely  based 
technology  rather  than  local  indicators  and  have,  to  a  great 
degree,  removed  the  "art”  from  single-stat ion  forecasting. 

However,  large-scale  conflict  or  operational  considera¬ 
tions  could  deny  all  forecasting  inputs  available  to  the 
forecaster  except  local  indications.  Reliance  on  single¬ 
station  forecasting  for  naval  operations  during  times  of 
crisis  or  war  was  first  addressed  by  Oliver  and  Oliver 
(1945)  : 


During  the  last  few  years,  wartime  conditions  have  made 
it  necessary  for  isolated  combat  units  to  issue  fore¬ 
casts  in  regions  where  no  network  of  meteorological 
stations  could  be  available.  Frequently  the  data  from 
several  stations  or  from  reconnaissance  planes  are 
available,  but  in  sane  regions  the  forecaster  must  rely 
only  on  surface  and  upper-ai>-  observations  made  at  his 
own  station.  This  is  particularly  true  in  the  case  of 
ships  at  sea.  Hence,  it  is  important  to  develop  profi¬ 
ciency  in  extracting  information  from  limited  aerolo- 
gical  data. 


While  the  above  quote  pertains  to  atmospheric  forecasts,  the 
same  could  be  true  for  ocean  predictions. 

In  modern  warfare,  isolated  ships  or  battle  groups  far 
from  the  main  force,  operating  in  an  emission  control 
(EHCON)  posture,  must  be  able  to  predict  their  operating 
environment  without  shore  support.  Such  predictions  are 
important  for  the  effective  utilization  of  weapon  and  sensor 
systems  and  to  deny  the  effective  use  of  the  environment  to 
the  opposing  force.  Affected  systems  include:  navigation 
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is  the  basis  for  an  entrainment  hypothesis.  A  closed  system 
of  equations  is  obtained  by  using  the  bulk  buoyancy  and 
momentum  equations  with  the  mean  turbulent  field  modeling  of 
the  vertically  integrated  equations  for  the  individual  IKE 
componen  ts. 

To  better  define  the  mixing  process,  separate  vertical 
and  horizontal  equations  for  THE  are  used.  Buoyancy  flux 
and  shear  production  provide  energy  for  vertical  mixing, 
with  buoyancy  flux  being  a  somewhat  more  efficient  source  of 
energy  for  mixing  due  to  its  direct  contribution  to  the 
vertical  component  of  TKE.  The  buoyancy  equation  is  derived 
from  the  heat  and  salt  equations  using  a  linearized  equation 
of  state: 


P  »  Por  1  -  «  (9  -  ao )  -  0  (S  -  S0)  ]  (2.  15) 

Buoyancy  is  given  by: 

b  =  g  (po  "  P  )  /P0  (2.16) 

where : 

Q  =  temperature 
S  =  salinity 
g  =  gravity 
P  =  density 

«  =  thermal  expansion  coefficient 

0  =  density  coefficient  for  salt 

The  tilde  represents  the  total  instantaneous  value  and  the 
subscript  "o"  denotes  an  arbitrary,  but  representative, 
constant  value.  In  the  short  term,  salinity  has  little 
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where  a  is  Stefan’s  constant  (4.61  X  10'11)  and  €c  is 
obtained  from  equation  2.10  .  The  net  long  wave  flux  at  the 
surface,  Flonfl,  becomes: 


r  a£T.4  ~  €c T4  -  n  -€c)isl 


(2.  13) 


where  T  is  the  average  temperature  of  the  cloud. 

For  the  cloud  free  case,  the  net  fluxes  are  calculated 
at  the  top  of  the  mixed  layer  (Z  =  h)  and  at  the  surface  (Z 
=  0)  by  integrating  the  flux  eaissivity  profile  (Fleagle  and 
Businger,  1980),  which  is  a  function  of  the  water  vapor  and 
temperature  profiles.  The  net  long  wave  flux  at  the  surface 
for  the  clear  sky  case  is: 


Flong  Fu  Fd 


(2.  14) 


where  Fu  and  Fd  are  the  upward  and  downward  radiative  fluxes 
respectively. 


C.  OCEABIC  BOUHDAHY  LAYEH  (OBI)  MODEL 

The  OBI  model  is  a  one-dimensional,  second  order  turbu¬ 
lence  closure,  vertically  integrated  (bulk)  model  of  the 
upper  ocean  surface  turbulent  boundary  layer  or  mixed  layer 
developed  by  Garwood  (1977)  .  It  uses  the  continuity  equa¬ 
tion  for  an  incompressible  fluid,  the  first  law  of  thermody¬ 
namics  (heat  equation),  the  conservation  of  salt  equation, 
an  analytical  equation  of  state,  the  Navier-Stokes  equation 
of  motion  with  the  geostrophic  component  eliminated,  and  a 
two-component  vertically  integrated  TKE  budget. 

An  understanding  of  the  dynamics  of  the  entrainment 
process  is  necessary  to  predict  the  MLD  change  with  time. 
The  stable  water  mass  underneath  the  mixed  layer  is  destabi¬ 
lized  and  eroded  by  TKE  in  the  lixed  layer.  This  TKE  budget 
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aerosols).  Total  short  wave  radiative  flux  at  the  surface 
is  determined  by  computing  both  the  direct  and  diffuse  radi¬ 
ative  components.  Fairall,  et  al.  (1981)  gives  a  comprehen¬ 
sive  review  of  this  method.  Solar  zenith  angle  is  necessary 
in  order  to  compute  short  wave  radiative  flux  and  it  is 
determined  from  the  latitude#  julian  day  and  time  of  day  of 
the  initial  input.  The  fraction  of  short  wave  radiation 
from  the  sea  surface,  Ag,  which  is  an  important  short  wave 
ratiative  parameter,  is  prescribed  in  the  MABL  model  as  0.1. 

long  wave  radiative  flux  is  calculated  differently  for 
the  cloudy  and  cloud  free  cases.  In  the  cloudy  case,  the 
only  clouds  permitted  are  non-black  stratus.  Radiative  flux 
is  a  function  of  cloud  emissivity,  €c ,  which  depends  on  the 
total  liquid  water  content,  W,  cf  the  cloud.  W  profiles  are 
nearly  linear  with  height  above  the  LCL  (Davidson,  et  al. , 
1984)  and,  along  with  €c  can  be  described  by: 

W  =  0.5  pa(h  -  Zc)g,h  (2.9) 

€c=  1  -  exp(-aW)  (2.  10) 

where  pa  is  the  density  of  air  (1.25  X  10"3  gm/cm3)  ,  h  is 
the  height  of  the  mixed  layer  (cloud  top),  Zc  is  the  LCL 
(cloud  bottom),  a  =  0.158  (Slingo,  et  al.,  1982),  and  q|h  is 
the  liquid  water  content  at  the  cloud  top. 

Osing  the  Stef an-Eoltzmann  law,  the  net  long  wave  cloud 
top  radiation  flux,  lnh,  can  be  calculated  from  the  cloud 
top  temperature,  Th  ,  and  the  effective  radiative  sky 
temperature,  TJk  ,  cloud  bottom  net  radiative  flux,  Lnc,  is 
calculated  in  the  same  way  using  the  cloud  bottom  tempera¬ 
ture,  T  c  ,  and  the  SST,  T#.  These  fluxes  are  given  by: 

Ln„=€c<*(V  *  W)  {2*11} 

L„c  =  *  Tc4)  (2-  12> 
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T*  = 

C#^(0O-  9) 

(temperature  parameter) 

(2.4) 

a*  * 

cfl*(g0-  a) 

(moisture  parameter) 

(2.5) 

These 

fluxes  are  given 

by: 

u  *  w  ’ 

=  u** 

(momentum) 

(2.6) 

T*  w  * 

=  U*T* 

(sensible  heat) 

(2.7) 

g  •  w  * 

=  u*g* 

(latent  heat) 

(2.8) 

where  Cd  and  are  the  stability  dependent  wind  and  temper¬ 
ature  drag  coefficients,  9  is  the  potential  temperature  and 
g  is  the  specific  humidity.  The  subscript  "o'1  denotes  the 
surface  value. 

Because  radiation  is  a  primary  factor  in  the  OBI  devel¬ 
opment,  considerable  effort  has  gone  into  calculating  the 
radiation  budget  in  the  MAB1.  Uncertainty  in  background 
aerosols,  distribution  and  concentration  of  atmospheric 
absorbing  gases,  and  cloud  droplet  size  spectra  are  sources 
of  error  in  the  MABL  radiation  calculations.  The  short  and 
long  wave  fluxes  are  computed  separately.  Short  wave  radia¬ 
tive  flux  is  calculated  using  the  delta-Eddington  method 
(Joseph,  et  al.,  1976)  which  accounts  for  heating  of  the 
mixed  layer  by  solar  radiation.  Incident  flux  at  the  top  of 
the  mixed  layer  is  obtained  fro*  the  flux  at  the  top  of  the 
atmosphere  and  the  average  transmittance  in  each  of  15  bands 
covering  the  spectrum  from  0.2  um  to  1.7um.  Short  wave 
extinction  is  due  to  both  scattering  and  absorption.  Water 
vapor  in  the  atmosphere  is  the  primary  absorbing  constit¬ 
uent.  Scattered  radiation  for  is  a  second  short  wave  radia¬ 
tive  component  (diffuse  solar  radiation)  which,  in  the  mixed 
layer,  is  due  to  atmospheric  particles  (cloud  droplets  and 
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and  the  integration  of  the  moisture  budget  equation  (Q) 
method.  Gleason  (1982)  found  that  the  Q-method  has  the  most 
merit  and  was  used  in  this  thesis.  However,  this  method 
assumes  that  no  advection  will  take  place. 

Integrated  rate  equations  by  Tennekes  and  Dreidonks 
(1981)  are  used  to  predict  the  time  rate  of  change  of  the 
conservative  quantities  and  their  respective  jumps  at  the 
inversion.  The  vertically  integrated  equations  fcr  a 
conserved  property,  X,  are: 

h  (DX/Dt)  =  (w«  x*)  -(v7!7)  ♦  source  (2.1) 

o  n 

h(DaX/Dt)  =  hrx(3h/3t)  -  (v’T"')  +  (w'  x  * »  -  source  (2.2) 

o  n 

where : 

source  =  -  (Fn  -Fn  o  )/pCp  for  X  =  temperature 

source  =  0  for  X  =  humidity 

and  Gamma  ( r )  is  the  conservative  property  lapse  rate  above 
the  inversion  and  F„  is  the  net  radiative  flux.  The 
subscripts  "h"  and  Mo"  refer  to  inversion  height  and  surface 
values  respectively. 

The  entrainment  velocity  parameterization  is  based  on 
the  formulation  by  Stage  and  Bu&inger  (1981)  and  is  used  to 
close  the  system  of  equations  and  determine  the  time  evolu¬ 
tion  of  the  inversion  height.  The  closure  assumpton  is  that 
the  dissipation  rate  of  the  TKE  is  a  fixed  fraction  ( 1  —  A)  of 
the  production  rate.  "4"  is  the  entrainment  coefficient  and 
is  taken  to  be  0.2. 

Bulk  aerodynamic  formulas  described  by  Davidson,  et  al. 
(1984),  are  used  to  determine  the  surface  fluxes  of 
momentum,  sensible  heat  and  latent  heat.  They  are: 

u*  *  Cd^01o  (friction  velocity)  (2.3) 
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Fig.  2.4  Schematic  of  I spat.  Prescription  and 
Computing  Steps  in  UBL  Prediction. 


The  SABI  nodel  is  very  sensitive  to  the  prescribed  large 
scale  subsidence  which  noderates  the  depth  of  the  aoist 
marine  layer  and  hence,  cloud  formation.  Hethods  which  can 
be  used  to  conpute  the  subsidence  from  single  station  obser¬ 
vations  are:  the  kinenatic  aethod;  the  adiabatic  method; 
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B.  BAB I IE  ATH0SPHEB3C  BOON DA HI  LATER  (NABL)  MODEL 

Existing  integrated  MABL  models  are  zero-order,  two 
layer,  mixed  layer  models  consisting  of  a  well  mixed  turbu¬ 
lent  boundary  layer  underneath  a  relatively  non- turbulent 
free  atmosphere.  A  transitior  zone  (inversion)  of  zero 
thickness  (hence  the  term  zerc-order)  separates  the  two 
layers.  Because  it  has  no  thickness,  a  jump  occurs  at  the 
inversion  of  the  conservative  parameters.  The  model  was 
described  by  Stage  and  Businger  (1981),  who  formulated 
entrainment  energetics,  and  modified  by  Davidson,  et  al. 
(1984),  who  described  bulk  aerodynamic  formulae  for  the 
surface  layer. 

The  following  inputs  are  required  by  the  model: 

(1)  an  initial  atmospheric  sounding; 

(2)  the  mean  winds  at  a  level  within  the  surface  layer; 

(3)  the  surface  temperature;  and 

(4)  the  subsidence  rate. 

As  formulated,  ten  winds  can  be  input  over  the  forecast 
period.  The  sea-surface  temperature  (SST)  remains  unchanged 
and  the  surface  current  is  assumed  to  be  zero.  SST  and  the 
surface  current  change  in  the  coupled  version  due  to  the 
interaction  with  the  OBL  model.  Inversion  height,  mixed 
layer  values  of  temperature,  humidity  and  wind  are  predicted 
at  30-minute  intervals  in  both  the  coupled  and  uncoupled 
versions.  The  prediction  steps  are  shown  in  Fig.  2.4. 
Procedures  are  the  same  for  clear  and  cloudy  cases  except 
for  entrainment  computation  and  estimation  of  cloud  top 
cooling.  Formation  of  clouds  cr  fog,  cloud  top  cooling  and 
associated  entrainment  are  important  in  the  physical 
processes  in  the  NABL. 
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Fig.  2.1  Idealized  Atmospheric  and  Oceanic 
Boundary  layer  Tea p era tare  profile . 


BID  vill  shallow.  Fig.  2.2  shews  the  iaportance  of  mechan¬ 
ical  sizing  to  the  CBL  and  its  interaction  with  buoyancy 
fluzes.  A  primary  objective  of  this  thesis  will  be  to 
determine  the  effects  of  wind  stress  on  the  HID. 

Several  cause  and  effect  relationships  are  evident. 
First,  clouds  can  be  caused  by  changes  in  the  ocean  surface 
temperatures  which,  in  turn,  affect  the  radiation  budget  of 
the  OBL  which  results  in  surface  temperature  changes. 
Secondly,  the  stability  of  the  BABL  is  influenced  by  air- 
ocean  temperature  differences.  Again,  cloud  formation  plays 
an  important  role.  Such  relationships  make  a  coupled 
prediction  approach  necessary.  This  coupled  model  uses 
local  oceanic  descriptions  from  a  bulk  oceanic  sized  layer 
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extend  above  the  lifting  condensation  level  (LCL) ,  clouds  or 
fog  will  form  within  the  layer  and  greatly  affect  the  OBL. 

In  the  ocean,  the  OBL  is  defined  as  extending  from  the 
surface  to  the  top  of  the  seasonal  thermocline.  Within  the 
mixed  layer,  nearly  homogeneous  profiles  of  temperature  and 
density  are  evident.  Rhile  density  is  a  function  of  both 
salinity  and  temperature,  over  short  time  scales,  tempera¬ 
ture  has  a  far  greater  effect  (Miller,  1976).  At  the  bottom 
of  the  mixed  layer,  large  changes  in  temperature  and 
salinity  exist  with  increasing  depth.  This  is  the  top  of 
the  thermocline  region  and  serves  as  a  transition  from  the 
mixed  layer  to  the  usually  dynamically  stable  interior 
ocean.  Idealized  model  profiles  for  both  the  MABL  and  OBL 
are  depicted  in  Fig.  2.1. 

Radiation  is  another  source  of  energy  in  both  mixed 
layers.  While  most  solar  energy  will  penetrate  the  MABL, 
this  is  not  true  of  the  OBL.  Long  wave  radiation  is 
absorbed  and  emitted  by  the  first  few  millimeters  of  sea 
water.  Therefore  downward  turbulent  heat  flux  is  as  impor¬ 
tant  as  the  upward  flux.  Surface  heat  flux  in  the  OBL  is  a 
function  of  sensible  and  latent  heat  as  well  as  the  long  and 
short  wave  solar  radiation.  This  radiation  is  highly  depen¬ 
dent  on  cloud  formation  in  the  ABL  so  proper  ABL  modelling 
is  essential  for  realistic  ocean  NLD  predictions. 

While  buoyancy  forces  create  turbulence  in  both  boundary 
layers,  relatively  more  turbulence  in  the  OBL  is  produced  by 
wind  stress  (Davidscn  and  Garwcod,  1984) .  The  MLD  will 
deepen  if  the  wind  generated  downward  flux  of  momentum  and 
forced  convective  mixing  (of  surface  water  that  has  been 
cooled  by  radiative  heat  loss  or  evaporation)  provide  suffi¬ 
cient  TK2  to  erode  the  thermocline  and  entrain  cooler  water 
downward.  If,  on  the  other  hard,  the  near  surface  absorp¬ 
tion  of  solar  radiation  produces  a  downward  buoyancy  flux 
which  dampens  turbulence  and  provides  for  net  warming,  the 
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II.  agfifii  BACKGB001P  ANE  ANALYSIS  PBOCEDOBES 


A.  BOOIOABY  LAYEB  DESCBIPTION 

Thermodynamically  and  dynaxically  forced  marine  atmos¬ 
pheric  (MABL)  and  oceanic  (OBI)  boundary  layers  must  be 
considered  in  coupled  model  description.  The  air-sea  inter¬ 
face  has  adjacent  oceanic  and  atmospheric  turbulent  mixed 
layers  which  effectively  insulate  the  quasi-geostrophic 
oceanic  and  atmospheric  flows  above  the  top  of  the  MABL  and 
below  the  bottom  of  the  OBL.  The  primary  source  of  the 
turbulence  is  the  velocity  (shear)  and  buoyancy  (density) 
gradients  created  by  the  exchange  of  mass,  energy  and 
momentum  across  the  air-sea  interface.  Such  turbulence 
generates  turbulent  kinetic  energy  (TKE)  and  vigorous 
vertical  mixing  (Davidson  and  Garwood,  1984). 

Under  undisturbed  conditions,  the  MABL  consists  of  a 
moist,  well-mixed  layer  extending  from  the  surface  to  the 
capping  inversion.  At  the  inversion,  jumps  in  temperature 
and  humidity  occur  with  a  rapid  increase  in  temperature  and 
a  corresponding  decrease  in  humidity.  Within  the  MABL, 
equivalent  potential  temperature  and  specific  humidity  are 
conservative  quantities. 

A  consequence  of  the  large  velocity  fluctuations  and 
mixing  in  the  MABL  is  that  energetic  eddies  (extending  from 
the  surface  to  the  inversion)  entrain  warm,  dry  air  and 
bring  momentum  into  the  mixed  layer.  This  results  in  the 
upward  growth  of  the  layer.  Subsidence  arises  from  large- 
scale  atmospheric  forcing  and  tends  to  limit  the  depth  of 
the  mixed  layer.  Therefore,  the  change  in  the  inversion 
height  with  respect  to  time  is  a  function  of  the  entrainment 
rate  and  subsidence.  If  entrainment  causes  the  MABL  to 
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tested  and  evaluated  with  shipboard  applications  in  mind. 
Accurate  predictions  of  ocean  mix-^d  layer  depths  (NLD)  , 
thermocline  jumps  and  gradients  are  essential  in  optimizing 
anti-submarine  warfare  (ASH)  systems  and  employing  acoustic 
countermeasures.  These  parameters  are  dependent  on  the  heat 
fluxes  and  wind  stress  at  the  ocean  surface.  This  thesis 
will  examine  how  the  single-station  coupled  model 
predicticns  of  the  surface  fluxes  and  mixed  layer  depth 
compare  with  observations. 


radars,  communications  equipment,  over-the-horizon  (OTH) 
radars,  air  and  surface  search  radars,  missile  dnd  gun  fire 
control  radars,  laser  and  infra-red  (IH)  guided  missiles, 
hull  mounted  sonars,  towed  arrays,  and  variable  depth  sonars 
(VDS)  . 


Pig.  1.1  Components  of  the  Coupled  Atmospheric 
and  Oceanic  Boundary  Layer  flodel. 


A  coupled  atmospheric  and  oceanic  model  (O'Loughlin, 
1982)  has  been  developed  to  predict  ocean  and  atmospheric 
parameters  from  local  indications.  The  components  of  this 
coupled  model  are  depicted  in  Pig.  1.1.  This  model  is  being 
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effect  on  the  upper  ocean  density  profile  except  at  higher 
latitiudes  (Miller,  1976) .  This  leaves  temperature  as  the 
dominant  density  factor.  However,  by  using  buoyancy  (b) 
instead  of  only  temperature  (9)  allows  the  model  to  be 
applied  in  cases  where  evaporation  and/or  precipitation 
contribute  significantly  to  the  surface  buoyancy  flux. 

Model  initialization  requires: 

(1)  mixed  layer  temperature  and  salinity  profiles; 

(2)  wind-driven  horizontal  current  profiles. 

The  initial  salinity  profile  is  made  isohaline  if  unknown 
and  the  initial  currents  are  assumed  to  be  zero  if  unknown. 
Lack  of  initial  information  about  current  and  salinity 
profiles  is  not  a  serious  deficiency  because  the  model  will 
evolve  reasonable  transient  profiles  of  salinity  and 
momentum  after  only  a  single  diurnal  cycle,  and  the  final 
results  are  not  sensitive  tc  the  initial  current  and 
salinity  profiles  (Davidson  and  Garwood,  1984).  It  is  also 
possible  to  prescribe  an  upper  ocean  internal  vertical 
velocity  (upwelling  and  downwelling)  if  it  is  known  to  be 
significant.  The  mixed  layer  depth,  h,  is  defined  as  the 
shallowest  depth  at  which  the  observed  density  value,  ot,  is 
0.02  ut  units  greater  than  the  observed  surface  density 
value. 

At  each  one  hour  model  time  step,  the  following  boundary 
conditions  are  required: 

•  net  upward  turbulent  heat  flux  at  the  water  surface 
(sensible  and  latent  heat  flux)  plus  back  radiation; 

•  incident  solar  radiation; 

•  the  fraction  of  short  wave  radiation  absorbed  in  the  top 
one  meter  of  the  ocean; 

•  surface  wind  speed  and  direction; 
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•  clcud  cover; 


•  temperature  (SST,  dry  bulb  and  dew-point  air  tempera¬ 
ture)  ;  and 

•  precipitation  (P). 

Vertical  velocity  at  the  botton  of  the  mixed  layer  may  be 
prescribed  hourly.  Other  physical  and  model  constants  that 
need  to  be  prescribed  are: 

•  the  extinction  coefficient  for  short  wave  radiation 
absorption; 

•  the  Coriolis  parameter;  . 

•  the  critical  Richardson  number  for  stability  adjustment 
below  the  mixed  layer; 

•  expansion  coefficient  for  temperature;  and 

•  density  coefficient  for  salt. 

The  model  forecast  is  not  particularly  sensitive  to  these 
constants,  but  reasonable  values  have  been  determined  by 
Gallacher,  et  al.  (1983),  and  they  need  not  be  readjusted 
for  geographic  and  seasonal  variability. 

By  use  of  bulk  aerodynamic  formulas,  the  turbulent 
fluxes  of  latent  heat,  Qe  ,  and  sensible  heat,  Qh  ,  can  be 
estimated  as  follows: 

Q.»  Cd  <.98Es  -  E.)010  (2.17) 

2h=  Cd(T,  -  Ta)01o  (2.18) 

Ret  back  radiation  is  estimated  from  the  empirical  equation 
(Husby  and  Seckel,  1S78); 
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Q  b  =  1.14  X  10-^(273.  16  ♦  Ts)  ♦  (.39  -  .5E#V2)  (2.19) 

(1  -  .  6  C2) 

where: 

Es  =  saturated  vapor  pressure  of  the  marine  air  (0.98 
corrects  for  salt  defects) 

Ea =  vapor  pressure  of  air  based  on  dew  point  temperature 

Ta=  air  temperature 

T, =  sea  surface  temperature 

C  =  fractional  cloud  cover 

Cd=  drag  coefficient 

The  upward  heat  flux,  Qu,  is  then  given  as: 

2U  =  2e  +  2h  +  Qb  (2.20) 

and  the  solar  radiation,  Qs,  by: 

Q,  “  (1  -  a«b  )  (1  -  . 66C3 )  Q0  (2.21) 

The  constants  "a"  and  "b"  are  adapted  from  Tabata  (1964)  and 
the  cubic  cloud  cover  correction  from  Laevastu  (1960). 
Clear  sky  radiation,  Q0,  is  given  by  Seckel  and  Beaudry 
(1973)  : 

(2.22) 

were  calculated  by 
predicted  in  the 


Q„  =  Aq  ♦  A,  cos  0  ♦  Bt  sin  0  +  A2cos  20 
♦  B2sin  20 

where  the  coefficients  (A0,  At  ,  etc.) 
harmonic  representation  of  the  values 


(2. 23) 


Smithsonian  Meteorological  Tables  (List  1958)  with: 

(J>  =  (2V365)  (t  -  21) 

where  "t"  is  the  julian  day  of  the  year. 

The  ocean  is  not  very  transparent  to  solar  radiation. 
In  the  open  ocean,  approximately  50%  is  absorbed  within  the 
first  meter  (this  fraction  is  called  "RfM).  Absorption 
varies  from  region  to  region  and  is  highly  dependent  on  the 
concentration  of  absorbing  particles  such  as  phytoplankton, 
"Gelbstoff”  (yellow  substance) ,  and  suspended  particulate 
matter.  Coastal  regions  absorb  more  radiation  because  of 
increased  amounts  of  suspended  particulates.  Very  little  of 
this  absorbed  radiation  penetrates  below  the  mixed  layer 
because  of  the  low  thermal  conductivity  of  the  underlying 
stable  thermocline.  Thus  most  cf  this  energy  is  transferred 
upward  out  of  the  ocean  and  back  into  the  atmosphere.  What 
short  wave  radiation  remains  penetrates  the  mixed  layer  and 
is  attenuated  in  an  exponential  fashion  depending  on  water 
turbidity.  Net  heat  flux  at  the  surface  is  given  as: 

Qn.t  *  Cu  ♦  (Rf)Q,  ~  Q,  (2.24) 

Surface  buoyancy  (heat  and  salt)  and  momentum  fluxes  can 
be  computed  by  use  of  the  foregoing  eguations.  Mixed  layer 


turbulent 

fluxes  are 

temperature,  salinity, 

given  by: 

velocity  and  buoyancy 

( T '  w ' )  = 

-Qn.t  /PCp 

(2.25) 

(SV)  = 

(P  -  E )S0 

(2.  26) 

(U’¥»)  = 

U** 

(2.27) 

(irnr*)  = 

g[«  (f^T7)  -  /3  (sTT7)  ] 

(2.  28) 
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where  the  subscript  no"  refers  to  the  surface  value.  The 
friction  velocity  in  air,  u*,  is  calculated  by: 

rs  =  P.Cd°io2  (2- 29) 

U*  =  (%/p.)  V*  (2.30) 

where  r  is  the  surface  stress  (dynes/cm*) .  The  fluxes  of 
momentum,  radiation,  latent  and  sensible  heat  at  the  sea 
surface  deteraine  shallowing  (retreat)  or  deepening  of  the 
mixed  layer  by  entrainment.  If  there  is  a  positive  buoyancy 
flux  (Qn#t<0  and  E>P)  the  MLD  will  deepen.  This  will  be  the 
case  at  night  or  whenever  long  wave  radiative  cooling  plus 
the  upward  turbulent  fluxes  of  heat  and  moisture  exceed 
solar  radiation.  Negative  buoyancy  flux  will  result  in  the 
shallowing  of  the  MID  and  is  caused  by  the  domination  of 
daytiae  solar  heating  (net  warning)  at  the  surface,  provided 
wind  stress  is  not  large. 

A  schematic  of  the  ocean  acdel  computations  is  shown  in 
fig.  2.5.  New  ocean  temperature,  salinity  and  wind-driven 
current  profiles  are  predicted  at  one-hour  intervals  and  are 
used  to  predict  the  new  MLD. 

D.  COUPLED  BO 01 DAB Y  LAYEB  MODE I 

Coupling  of  the  ataospheric  and  oceanic  models  was  first 
accomplished  by  O'Loughlin  (1982)  by  matching  momentum, 
sensible  heat,  latent  heat,  and  radiation  at  the  air-sea 
interface.  The  ataospheric  part  of  the  coupled  model 
computes  the  radiative,  heat  and  momentum  fluxes  which  are 
inputs  to  the  ocean  model.  Since  the  ocean  model  has  an 
hour  time  step,  and  the  ataospheric  model  has  a  half-hour 
tine  step,  the  ocean  model  was  inserted  as  a  subroutine  into 
the  ataospheric  model  and  is  called  every  other  atmospheric 
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time  step.  Before  the  atmospheric  flaxes  can  be  utilized  by 
the  ocean  model,  they  aust  be  converted  to  the  appropriate 
units  for  seawater. 

The  ocean  model  prediction  requires  the  vector  surface 
stress  for  the  ocean  turbulent  velocity  flux,  u*w*.  This 
requires  that  the  wind  field  be  broken  into  horizontal 
coaponec  ts: 


Oio„  =  ”  (si®  ®)  ®io  (2. 

°ioy  =  "<cos  9)°10  (2.  32) 

where  01Q  is  the  wind  speed  at  10  meters  and  01Ox  ,  u10y  are 
the  east-west  and  north-south  horizontal  components  respec¬ 
tively.  The  wind  direction,  9,  is  relative  to  true  north. 
Then  u*w2  is  given  by: 

r.  ■  P.u*2  *  P*u*w2  (2.33) 

the  density  of  seawater,  pw,  to  be  about  1  gm/cm3.  Finally, 
the  horizontal  components  of  the  turbulent  velocity  fluxes, 
u*wx 2  and  u*wy2,  are  computed  by: 


U*2  = 

C*°102  =  V°10x2  *  °10y  2) 

(2.  34) 

U*  2 

*1 

*  VlOx  2 

(2.35) 

U*Wy2 

*  Vioy  2 

(2.  36) 

SST  is 

passed  from  the  MABL  mcdel  back  to 

the  atmospheric 

model  for  use  in  the  next  time  step.  Fig.  2.6  shows  the 
interrelationships  between  the  atmospheric  and  oceanic 
models  in  the  coupled  model. 

The  coupling  process  has  the  potential  of  significantly 
enhancing  the  OBI  prediction  qualities  because  of  the 
improvement  in  the  boundary  conditions  inherent  in  a  coupled 
OBL-HABL  system.  longer  period  reasonable  forecasts  should 
be  obtained  due  to  the  feedback  between  the  adjacent  turbu¬ 
lent  boundary  layers  and  the  associated  thermodynamic  and 
dynamic  adjustments  in  each  layer.  Formation  of  stratus  or 
fog  in  the  HABL  model  is  particularly  important.  Accurately 
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prescribed  incoming  solar  radiation  is  inperitave  for  a 
useable  CBL  forecast.  This  can  only  be  done  by  use  of  an 
MABL  model  such  as  has  been  described  in  section  B  (Davidson 
and  Garwood,  1984). 


E.  Q*  ESTIMATION  FBOH  T0BBD1ENT  KINETIC  ENEBGY  (TKE) 
DISSIPATION 

Direct  measurement  of  u*  on  a  ship  is  rather  difficult. 
One  popular  technique  which  is  used  extensively  at  the  Naval 
Postgraduate  School  is  the  dissipation  method.  An  excellent 
review  of  this  method  as  well  as  the  bulk  method  is 
contained  in  Schacher,  et  al.  (1981).  The  dissipation 
method  involves  the  use  of  the  turbulent  kinetic  energy 
(TKE)  equation  which,  assuming  steady  conditions  and  the 
divergence  term  to  be  small  relative  to  the  other  terms, 
becomes; 

u*  =  {l€kz]/[<t>m(z/L)  -  (z/L)  ]}V*  (2.37) 

The  denominator  on  the  right  hard  side  of  equation  2.37  has 
teen  parameterized,  where  the  quantity  4)m(z/L)  -  (z/L)  is 

substituted  by  <J>e  following  Nyngaard  and  Cote  (1971): 

<J>€  »  (1  ♦  0.5  |z/L|  =¥**)#  for  Z/L  <  0  (2.38) 

4»€  =  (1  ♦  2.5  (z/I)  ^)  ^*  for  z/L  >  0  (2.39) 

Determination  of  the  dissipation  rate,  €,  necessary  for  this 
technique  requires  that  the  inertial  subrange  follows  a 
Kolmogorov  spectrum,  i.e. ,  where: 

Su  (k*)  =  «  e^*k*-V3  (2.40) 

Su(k*)  in  equation  2.40  is  the  spectral  density  of  the  hori¬ 
zontal  windspeed,  «  is  a  coefficient  with  a  magnitude  of 
0.52  (Champagne,  et  al. ,  1977),  and  k*  is  the  wavenumber. 

Equation  2.40  has  been  experimentally  verified.  On  the  R/V 
Acania  a  hot  film  sensor  was  used  to  obtain  a  spectrum  of 
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high  frequency  velocity  fluctuations  from  which  the  IKE 
dissipation  u*  was  calculated.  It  is  important  to  note  here 
that  poor  velocity  turbulence  spectra  are  obtained  during 
periods  cf  fog,  rain  cr  drizzle.  This  leads  to  the  dissipa¬ 
tion  method  being  invalid.  Errors  may  also  occur  if  the 
relative  wind  is  less  than  3  m/s  allowing  condensation  to 
form  on  the  hot  film. 

Applying  Taylor's  hypothesis,  i.e.,  k*  =  2*f/0,  equation 
2.40  may  be  rearranged  to  become: 

€  =  240  [Su(f)  ]0  (f/0)  0  (2.41) 


where  f  is  frequency  measured  in  Hz.  The  inertial  subrange 
spectral  densities  are  integrated  between  two  wavenumbers,  k^ 
and  kh  respectively,  in  equation  2.40  following  Khalsa  and 
Bussinger  (1977) : 


°Ak2 


(k)  dk  =  (3/2) 


«€0(k/0  -  k"0) 


(2.42) 


Applying  Taylor's  hypothesis,  equation  2.42  now  becomes: 


orAfa  =  0.230  (€0)0  (f/0  -  f’0)  (2.43) 

Rearranging  equation  2.43  and  combining  with  equations  2.37, 
2.38  and  2.39,  with  a  von  Karman  constant  of  0.4,  we  obtain: 


u*  =  2.81  cru  (z/[U<J)  (z/1)  ]}  0  (f/0  -  f‘0)  0  (2.44) 

A  f  c 

Substituting  f  =  50Hz,  f  =  5Hz,  the  measurement  height,  z  = 
20.0m,  equation  2.44  reduces  to  its  operational  form  as 
follows: 


u*  =  3.69  *v[02OiA(z/L>  ]"* 
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(2.45) 


/ 


Because  z/L  is  a  function  of  the  friction  velocity,  u* 
according  to  equations  2.38  and  2.39,  the  solution  to  egua 
tion  2.45  involves  an  iterative  process  (Large,  1979). 


III.  SYHOPTIC  SITUATION 


A.  SIZED  LAYEH  DY1ABICS  EX  PEE  I  BENT  (MIIDEI) 

The  Mixed  layer  Dynamics  Experiment  (MILDEX)  was 
conducted  during  the  period  of  24  October  to  10  November, 
1983  near  34°N  and  126°W.  MILDEX  was  a  multi-group  and 
aulti-platf ora  experiment  designed  to: 

•  provide  magnitudes  of  air-sea  energy  exchange  rates  for 
use  as  boundary  conditions  in  mixed  layer  modelling. 

•  evaluate  the  drag  coefficient  from  TXE  dissipation  meas¬ 
urements  and  its  dependence  on  swell  amplitude  and 
direction. 

•  to  provide  a  time  series,  from  radiosonde  measurements, 
of  MAB1  structure  for  model  verification. 

•  provide  an  intercom pari son  between  radiative  transfer 

algorithms  developed  at  Scripps  Institute  of 

Oceanography  and  aeasured  data. 

•  evaluate  the  utility  of  SODAR  as  a  ship-torne 

instrument. 

•  evaluate  radiative  transfer  models  of  the  atmosphere 
with  cloudiness  as  a  unique  parameter. 

Meteorology  data  were  collected  on  the  R/V  Acania  as  well  as 
two  other  platforms  during  MILDEX.  R/V  Acania  data  collec¬ 
tion  consisted  of  measurements  of  windspeed,  temperature, 
humidity,  radiation  and  atmospheric  pressure  in  the  surface 
layer,  and  sea  surface  temperature  (measured  by  thermistor)  , 
and  radiosonde  derived  profiles  of  atmospheric  windspeeds 
and  directions,  teiperatures  and  humidity.  Radiosonde 
launches  were  coordinated  with  satellite  pass  times. 
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Fig.  3.1  depicts  the  locations  on  the  bow  of  the  R/V 
Acania  of  mean  and  turbulent  aeasurement  sensors.  There 
were  two  levels  of  these  measurement  sensors:  bow  mast  at  5 
meters  and  main  mast  at  20  meters.  In  addition,  there  were 
hourly  observations  of  sea  state  and  cloudiness. 
Conductivity,  temperature  and  depth  (CTD)  casts  were  also 
made  from  the  B/V  Acania.  Two  ether  vessels,  the  R/V  Wecoma 
and  R/P  FLIP  collected  meteorological  and  oceanic  data  near 
the  same  site  as  B/V  Acania.  During  the  periods  of  interest 
in  this  thesis,  the  R/P  FLIP  was  always  within  two  kilo¬ 
meters  of  the  B/V  Acania.  The  meteorological  data  has  been 
compiled  in  the  MILDEX  report  (Geernaert,  et  al. ,  1983)  . 


MET.  TOWERS 

MEAN  AND 
TURBULENCE 
SENSORS 


Fig.  3.1  Schematic  of  the  B/V  Acania 
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Time  series  of  wind  speeds,  directions,  SST,  air  temper¬ 
ature,  and  surface  pressure  in  Fig.  3.2  show  seme  of  the 
variability  experienced  during  aiLDEX.  Frontal  passages  are 
evident  cn  29  October,  31  Octcber,  and  06  November.  For 
most  of  the  experiment,  the  SST  was  greater  than  the  air 
temperature  (except  cn  29  October  and  09  November) ,  so  the 
surface  layer  was  unstable  or  near  neutral.  Only  for  6-8 
hour  periods  on  29  October  and  09  November  was  the  surface 
layer  stable. 

Three  periods  (cases)  were  considered  for  purposes  of 
this  study.  These  cases  were  chosen  on  the  basis  of  strong 
wind  forcing  of  the  OBL  and  the  availability  of  reliable 
radiosonde  data  for  model  initialization.  Fig.  3.3  shows 
the  movement  of  the  R/V  Acania  for  all  three  cases.  The  R/V 
Acania  was  located  near  34. ON,  125. 4N  during  all  periods. 
All  times  are  Pacific  Daylight  Time  unless  otherwise  speci¬ 
fied.  NABL  turbulent  kinetic  energy  dissipation  rate 
derived  u*'s  were  obtained  from  velocity  fluctuation  meas¬ 
urements  from  the  R/V  Acania  main  mast  (at  20  m)  and  a  bow 
mast  (at  5  m)  instruments.  MID  data  were  obtained  from  CTD 
casts  from  the  R/V  Acania  for  Case  I  and  by  the  R/P  FLIP  for 
Cases  II  and  III. 

Conditions  associated  with  these  cases  was  dominated  by 
increasing  clouds,  swell  and  winds  with  intermittent  to 
heavy  rainfall.  Heavy  rain  showers  occurred  on  29  October 
when  the  first  of  several  fronts  passed  through  the  area. 
Since  the  presence  of  fog,  rain  or  drizzle  affects  the 
validity  of  the  hot  film  signals  used  in  the  dissipation 
technique,  care  was  taken  to  exclude  affected  results  during 
the  periods. 

Ocean  conditions  were  also  variable.  MILDEX  experiment 
participants  reported  large  scale  advection  and  upwelling 
occurred  during  the  course  of  the  experiment.  Significant 
internal  wave  activity  was  also  evident  and  is  being 
analyzed  by  other  investigators. 


Speed  (M/S) 


TABLE  II 

Case  One  Averaged 

Binds 

Time  (P DT) 

S£eed  (n^s) 

Direction 

0900 

8.  3 

181.0 

1130 

8.9 

178.0 

1400 

8.  2 

203.  0 

1630 

8.6 

196.  0 

1900 

4.6 

263.  0 

2130 

5.4 

221.0 

0000 

5.7 

253.  0 

0230 

4.5 

253.  0 

0500 

4.2 

231.0 

0730 

4.9 

295.0 
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veraged  Hinds. 


B.  CASE  I  {0900,  31  OCTOBER  -  0900,  01  NOVEMBER) 

This  first  case  resulted  in  the  poorest  comparison  of 
predicted  and  observed  values  of  any  of  the  three  cases. 
Initialization  data  is  summarized  in  Table  I.  Table  II 
lists  the  2.5  hour  averaged  winds  that  were  used  for  subseq¬ 
uent  runs.  In  Fig.  4.1,  a  comparison  of  the  two  wind  inputs 
is  shown.  Fig.  4.2  shows  a  comparison  of  u*'s.  A  diffi¬ 
culty  in  comparing  model  and  instrument  derived  values  for 
this  case  was  instrument  performance  and  equipment  break¬ 
down.  Both  of  these  problems  were  encountered. 


Rain,  drizzle,  and  low  relative  wind  speeds  affected  the 
accuracy  of  both  the  bow  and  main  mast  dissipation  derived 
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IV.  DATA  &NI  RESOLTS 


A.  APPROACH  TO  FOBCIBG  0*  INTERPRETATION 

For  each  of  the  three  cases,  the  u*  associated  with  the 
model  predictions  was  compared  with  the  TKE  dissipation  rate 
{€)  derived  u*  and  the  model  MLD  was  compared  with  the 
observed  MID.  As  described  previously,  the  epsilon  (e) 
estimates  necessary  fcr  u*  calculations  by  this  method  were 
obtained  from  hot  film  measurements  on  the  bow  (5  m)  and 
main  masts  (20  m)  of  the  R/V  Acania  (Fig.  3.  1)  . 

The  MA3L  model  is  configured  for  ten  input  winds,  so  the 
ten  greatest  and  least  winds  were  chosen  from  each  24-hour 
period.  These  will  be  referred  to  in  the  discussion  as 
"peak”  winds.  Linear  interpolation  was  used  for  data  points 
between  these  ten  points.  However,  a  forecaster  at  sea 
would  not  predict  high  and  low  winds,  but  would  most  likely 
forecast  average  winds  over  the  period.  Therefore,  another 
model  run  utilizing  2.5  hour  wind  averaging  was  input  over 
the  same  periods  to  see  if  the  use  of  averaged  winds  would 
have  any  deleterious  effect  on  the  MLD  computations. 

The  purpose  of  these  comparisons  is  to  evaluate  atmos¬ 
pheric  forcing  (u*)  in  relation  to  predicted  and  observed 
MLD  *  s .  There  are  three  u*'s  used  in  this  study:  those 
estimated  by  bulk  formulations  using  peak  winds,  by  bulk 
formulations  using  averaged  winds,  and  by  dissipation  esti¬ 
mates.  In  the  following  discussion,  the  first  two  sets  will 
be  referred  to  as  model  u*  (peak  or  averaged)  and  the  last 
set  as  dissipation  u*. 
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Fig.  3.7  GOES  SEST  Visual  2215  GMT  NOVEMBER  03,  1983 


■HC  Surface  Analysis  1200  GMT  NOVEMBER  03,  1983 


GOES  SEST  Visual  2215  GMT  NOVEMBER  01,  1983 


WEST  visual  imagery  for  Case  III.  At  the  beginning  of  the 
period,  50%  cloud  cover,  consisting  mainly  of  stratocuaulus, 
was  evident.  This  changed  to  cirrus  and  finally  a  clear 
night-tiae  sky.  On  the  morning  of  5  Nov.  there  was  30*  to 
50S  cloud  cover  (mostly  cirrus  and  some  cumulus)  with 
clearing  in  the  afternoon.  Nc  precipitation  was  recorded 
for  this  period  which  was  good  for  hot  film  performance. 
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Swells  were  significant  and  fros  the  NW  and  ranged  from  7  to 
9  ft.  Fluctuations  in  air  temperature  never  exceeded  1.0°C. 
As  in  Case  I,  the  SSI  was  greater  than  the  air  temperature 
with  the  air-sea  temperature  difference  varying  from  -.6°to 
-1.6  °C.  In  this  case,  the  MAEL  surface  layer  was  unstable 
throughout  the  entire  period. 

Fig.  3.7  shows  GOES  WEST  visual  imagery  for  this  case. 
Cloud  cover  was  heavy  in  the  late  afternoon  and  again  at 
mid-morning  on  the  next  day.  Curing  the  night,  the  sky  was 
cloud-free.  Fog,  cumulus  and  stratus  clouds  were  evident  in 
the  afternoon.  Stratus  with  some  cirrus  was  present  the 
next  morning.  No  precipitation  was  recorded  during  the  Case 
II  period. 

D.  CASE  III  (1600,  Oh  NOVEflBEB  -  1600,  05  M07EHBEB) 

The  cold  front  approaching  during  Case  II  dissipated  at 
the  beginning  of  the  period  and  the  area  came  under  the 
influence  of  a  weak  high  pressure  south  of  the  R/V  Acania. 
Fig.  3.8  is  the  NMC  surface  analysis  for  Case  III.  Another 
low  pressure  system  is  to  the  NW  with  an  associated  cold 
front.  This  second  system  did  not  appear  to  influence  the 
weather  or  the  seas  in  this  case.  Winds  initially  were  from 
the  NNW  and  gradually  shifted  to  the  SSW  as  the  high  moved 
through  the  area.  Speeds  generally  decreased  from  about  7 
to  2  m/s,  and  air  temperature  fluctuated  no  more  than  .7°C. 
Surface  pressure  varied  by  2  mb  from  1022  to  1024  mb  while 
sea  swells  were  from  the  NW  at  8  to  10  ft.  As  in  the 
previous  two  cases,  the  5ST  was  greater  than  the  air  temper¬ 
ature  with  the  air-sea  temperature  difference  varying  from 
-1.  1  to  -2.1°C.  Again,  the  J1ABI  was  unstable  for  the  entire 
period. 

This  last  case  had  the  least  cloud  cover  due  to  the 
influence  of  the  high  pressure  system.  Fig.  3.9  shovs  GOES 
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was  located  in  the  northeast  Pacific  ocean,  and  a  high  domi- 
nated  the  Pacific  north  of  Hawaii.  The  entire  period  had 
nearly  1005E  cloud  cover  and  intermittent  rain  or  drizzle. 
Several  rain  squalls  were  noted  after  the  frontal  passage. 
Clouds  varied  from  altostratus  prior  to  the  frontal  passage 
to  stratus  for  the  rest  of  the  Case  I  period.  No  fog  was 
reported,  so  a  considerable  amount  of  hot  film  data  was 
available.  Fig.  3.5  is  a  copy  cf  GOES  WEST  visual  imagery. 

As  shown  in  Fig.  3.2,  winds  were  initially  from  the  SSW 
with  a  shift  to  the  west  after  the  frontal  passage  with 
speeds  varying  from  7  to  9  m/s  ahead  of  the  front  and  from  3 
to  5  m/s  behind  the  front.  There  was  a  steady  decrease  in 
surface  pressure  from  1018  mb  tc  1014  mb.  Swells  increased 
from  4  to  9  ft  after  the  frontal  passage  and  gradually 
decreased  to  5  ft  towards  the  end  of  the  observational 
period.  Air  temperature  dropped  by  1.5°C  after  the  frontal 
passage.  Throughout  the  period,  the  SST  was  greater  than 
the  air  temperature  with  the  greatest  difference  at  night  as 
expected.  Air-sea  temperature  differences  varied  from  -.2° 
to  -2.3°  C.  The  MABL  surface  layer  was  near  neutral  to 
unstable  during  the  period. 

C.  CASE  II  (1 600,  03  HOYEHBEB  -  1600,  04  NOTEHBEB) 

During  this  period,  a  stronger  cold  front  than  in  Case  I 
was  approaching.  Fig.  3.6  is  the  NMC  surface  analysis  for 
Case  II.  A  low  center  was  positioned  in  the  Gulf  of  Alaska, 
and  a  high  pressure  center  had  intensified  west  of  Baja 
California  and  south  of  the  Acania.  No  frontal  passage  was 
observed  for  this  case.  Binds  were  generally  f the  NNW 
shifting  to  the  west  in  the  middle  of  the  observational 
period  and  then  returned  to  the  NNH  at  the  end  while  wind 
speeds  fluctuated  between  3  and  7  m/s.  The  surface  pressure 
was  very  steady  changing  only  bj  2  mb  from  1021  to  1023  mb. 
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Pig.  3.3  Soweaent  of  H/f  Acania:  Case  I  (solid) , 

Case  II  (dash)  ,  Case  III  (dot)  . 

B.  CASE  I  (0900 ,  31  0CT0BEB  -  0900,  01  NOYENBEB) 

During  this  period,  the  H/V  Acania  experienced  weak  ccld 
frontal  passage  at  about  1700  or  31  October.  Fig.  3.4  shows 
the  National  Meteorological  Center  (NSC)  surface  analysis 
for  Case  I.  The  synoptic  situation  was  that  a  low  center 
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Fig.  4.2  Case  I  Friction  Velocities 

u*  values.  Post-processing  removal  of  the  bad  data  resulted 
in  significant  data  gaps.  For  the  bow  mast,  dissipation  u* 
values  are  missing  from  1616  tc  2315  and  from  0604  onward. 
Only  12  useable  observed  data  points  were  available  (as 
compared  to  48  for  the  model)  fcr  the  period.  The  field  log 
indicates  a  lack  of  confidence  in  main  mast  spectra  by  the 
observers  onboard  due  to  spurious  voltage  readings  and 
connector  problems,  so  data  from  it  were  not  included  at  all 
in  this  comparison. 

Fig.  4.2  shows  the  peak  wind  and  averaged  wind  model  u* 
and  the  dissipation  u*  values.  As  expected,  the  model  u* 
values  are  a  function  of  the  differences  in  the  input  winds. 
Fhere  the  averaged  winds  are  higher,  the  u*'s  are  higher. 
The  dissipation  u*  values  obtained  from  the  bow  mast  have 
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the  same  trend  but  are  lower  than  both  model  u**s  in  the 
pre-frontal  period  (before  17  0  0)  .  They  are  close  to  the 
"peaked"  wind  model  predictions  for  the  rest  of  the  time. 
All  of  the  u*  dissipation  values  are  lower  than  the  averaged 
wind  model  predictions. 

These  results  of  model  and  dissipation  u*  values  were 
not  completely  expected  for  the  frontal  regime.  Previous 
comparisons  of  bulk  and  dissipation  u**s  in  frontal  regions 
(Geernaert,  et  al. ,  1985)  indicated  that  the  dissipation 
u*'s  were  higher  than  bulk  values  prior  to  a  frontal  passage 
and  then  the  magnitudes  reversed  after  the  passage.  Here 
the  dissipation  u*’ s  remained  lower  than  the  bulk  value 
throughout.  Generally,  the  study  is  revealing  that  small 
pre-frontal  wind  speeds  (less  than  6  m/s)  produce  greater  u* 
deviations  (turbulent  vs  bulk) .  Additionally,  if  the  swell 
and  wind  oppose  each  other,  greater  u*  deviations  were  also 
evident.  These  were  minimized  if  the  wind  and  swell  were  at 
right  angles  to  each  other.  The  reasons  for  this  are  not 
entirely  clear.  However,  it  has  been  postulated  (Geenaert, 
et  al.,  1985)  that  the  direction  and  magnitude  of  the  swell 
relative  to  the  wind  (primarily  the  interaction  of  wind 
stress  with  short  wave  elements  on  long  wave  crests)  affects 
the  roughness  length,  zo,  which  in  turn  affects  u*  and  Cd . 

No  attempt  will  be  made  with  these  results  to  examine 
how  the  swell  interacted  with  the  wind  and  how  this  could 
have  modified  the  u*  values.  Perhaps  this  case  did  not  have 
the  same  kind  of  swell-wind  interaction  as  in  previous 
studies.  This  is  an  area  for  further  investigation. 

Fig.  4.3  compares  the  model  and  observed  HLD's.  The 
model  prediction  was  based  on  bulk  u*  values  shown  in  Fig. 
4.2.  In  this  case,  model  results  (solid  line)  early  in  the 
period  do  not  agree  with  the  observed  HID  (dashed  line) . 
While  both  curves  nearly  coincide  at  the  start,  the  model 
curve  shallows  much  more  quickly  although  it  reaches  the 
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Fig.  4.3  Case  I  Mixed  layer  Depth 

proper  depth.  Then  the  model  deepens  rapidly  and  stays 
significantly  deeper  for  the  rest  of  the  case.  An  MID  error 
of  up  to  10  meters  can  he  tactically  significant  for  naval 
forces  in  predicting  sonar  performance. 

The  early  shallowing  may  not  constitute  a  model 
prediction  error.  A  CTD  cast  would  give  the  thermal  MID 
while  the  model  predicts  the  turbulent  boundary  layer  depth. 
It  takes  a  finite  amount  of  time  for  the  temperature  field 
to  adjust  to  the  turbulent  field.  The  expression  for  the 
lag  is: 


t  ~  £  (aT)  u*3  ]/Q0  2 


(«.1) 


where : 

T  «  temperature 


t  =  lag  time 

Q0  =  surface  heating 

Thus  based  on  eguation  4.  1,  the  greater  the  u*,  the  greater 
the  lag  and  the  greater  the  heating  (Q0)  ,  the  less  the  lag. 

Osing  2.5  hour  averaged  winds  showed  no  improvement  in 
HID  prediction  (Fig.  4.3  dotted  line)  and  in  fact  gave  a 
poorer  prediction.  Uhile  the  two  wind  inputs  (Fig.  4.1)  are 
quite  different,  the  output  Mil's  are  quite  close  together 
except  at  1300.  This  is  a  pre-frontal  period  of  strong 
winds  and,  from  1200  to  1700,  the  averaged  wind  speeds  are 
higher  than  the  initial  winds.  Since  MID  is  a  function  of 
the  wind  speed  cu\)ed,  the  stronger  winds  would  give  deeper 
mixing.  This  is  again  true  for  the  period  from  2300  to  the 
end  of  the  run  where  the  averaged  winds  are  generally  higher 
than  the  initial  winds.  Overall,  one  can  conclude  that  for 
this  case  both  the  "peaked”  and  averaged  winds  give  poor  MID 
predictions  and  these  could  have  adverse  effects  for  naval 
operations.  However,  if  the  dissipation  u*  values  for  0900 
to  1700  (which  are  lower  than  the  bulk  values)  had  been  used 
by  the  model,  the  model  HID  would  have  been  closer  to  the 
observed. 

C.  CASS  II  (1600,  03  S07EHBER  -  1600,  04  N07EHBEB) 

The  coupled  model  MLD  prediction  agreed  well  with  the 
observed  HID  changes  in  this  case.  The  model  predicted 
unstable  ABL  conditions  throughout  the  period  which  agreed 
well  with  observations.  Table  III  contains  the  initializa¬ 
tion  data  for  this  case  and  the  2.5  hour  averaged  winds  are 
listed  in  Table  17.  Fig.  4.4  compares  the  two  input  winds 
for  this  case. 

This  was  a  period  with  unreliable  dissipation  u*  data. 
Fig.  4.5  is  a  plot  of  the  various  friction  velocities. 
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TABLE  HI 

Case  Two  Initialization 
Atmospheric  Values 


Pressure:  1022.4  mb 
Pot.  T(C):  17.1 
Cond 


Lift. 


Lvl:  476.71m 


Sea  Temp  (Cl:  18.78 
Spec.  Hum:  9.32 
Aavection:  0. 


Ocean  Values 


Sized  Layer  Depth:  12.0m  Temp.  Jump 
Below  Layer  Gradient:  -.0094 


(C):  -.2 


lije  (£j>£) 

.  1600 
1745 
2120 
2230 
2300 
0135 
0400 
0705 
1005 
1420 


Speed  (m/s) 

4.0 

5.6 

8.6 

4.7 

6.4 

2.4 
7.0 

hi 

4.6 


Direction 

040.0 

310.0 

036.0 

319.0 

355.0 

270.0 

281.0 

350.0 

008.0 

336.0 


Hhile  fog  or  rain  were  not  a  factor  in  dissipation  measure¬ 
ments,  equipment  failure  was  a  problem.  This  led  to  a 
reduction  of  dissipation  u*  data.  The  bow  mast  data  were 
not  available  from  1700  to  2300,  2330  to  0217,  and  0900  to 
1600  when  a  major  eguipment  repair  was  accomplished.  Bow 
dissipation  u*  data  from  12  thirty  minute  periods  were 
available.  Bain  mast  dissipation  measurements  were  not 
available  from  1600  tc  0752  and  from  0822  to  1340.  It  was 
also  absent  after  1410  and  only  data  from  4  thirty  minute 
periods  were  available. 

At  the  beginning  of  the  period,  the  bow  mast  dissipation 
u**s  are  close  to  the  bulk  derived  (model)  values.  However, 
during  the  middle  part  of  the  period,  the  bow  dissipation 
values  are  consistently  lover  than  the  bulk  (model)  although 
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Speed  (M/S) 


TABLE  IT 

Case  Two  Averaged 

Hinds 

i 

Time  {£DT ) 

S_geed  (m/s) 

Directi  on 

1600 

3.9 

040.0 

1830 

4.5 

345.0 

2100 

5.7 

349.0 

2330 

4.3 

308.0 

0200 

3.6 

274.0 

0430 

6.  1 

282.0 

0700 

4.6 

307.0 

0930 

6.6 

357.  0 

1200 

6.8 

348.0 

1430 

5.4 

341.0 

LEGEND 
INITIAL  WINDS 


1600  1300  2000  2200  0000  0200  0400  0600  0300  1000  1200  1400  1600 

Loco  Time  (Hours) 


f.  4.4  Coaparison  of  Case  II 
nitial  and  Averaged  Binds 


LEGEND 
MODEL  USTAR 


AVE  WIND  USTAR 

'  i" '  BOW  M  A  St  ‘OSTAR 
♦  MAIN  MAST  USTAR 


1600  1800  2000  2200  0000  0200  0+00  0600  0800  1000  1200  1+00  1600 

Local  Time  (Hours) 


Pig.  4.5  Case  II  friction  Velocity 

they  tend  to  have  the  same  time  variations.  The  two  bow  u* 
dissipation  values  near  0800  are  close  to  the  model  values. 

It  is  apparent  from  Fig.  4.6  that  the  model  predicted 
the  MID  well  enough  for  Naval  tactical  applications.  The 
observed  HLD  (dashed  line)  shots  some  differences  from  0200 
to  1000  which  are  probably  due  to  internal  wave  activity, 
advection  or  as  a  result  of  the  ship  moving  to  a  different 
area.  Observed  MLD  data  was  not  available  for  1100  to  1300. 
This  contributes  to  the  abrupt  MID  rise  at  1400.  The  model 
(solid  line)  predicted  the  MID  guite  well  especially  for  the 
deepening  part  of  the  period.  Better  agreement  between 
observed  and  model  shallowing  would  probably  be  evident  if 
MLD  data  would  have  been  available  for  1100  to  1300. 


Comparison  of  the  peaked  winds  with  the  2.5  hour  aver¬ 
aged  wind  HLD's  yields  some  interesting  results.  During  the 
deepening  phase,  the  averaged  winds  consistently  give  a 
shallower  MLD  (dotted  line)  than  the  initial  winds.  At  this 
tine  the  "peaked”  winds  are  consistently  stronger  (Fig.  4.4) 
than  the  averaged  winds  (except  at  0100)  and  would  produce 
more  vigorous  nixing.  However,  during  the  shallowing 
portion,  the  averaged  wind  prediction  is  significantly 
deeper  than  the  peaked  wind  prediction.  This  tine  the  aver¬ 
aged  winds  are  predominantly  stronger  than  the  peaked  winds. 
Again  the  predicted  MLD  results  are  expected  since  the  MLD 
is  proportional  to  u«.  if  the  dissipation  u*  values  for 
0200  -  0600  would  have  been  used  by  the  model,  the  predicted 
MLD  would  have  been  shallower  fcr  that  period  and  closer  to 
the  observed. 


D.  CASE  III  (1600,  04  N07EMBEB  -  1600,  05  H07EMBBB) 


As  in  the  previous  case,  the  coupled  model  HID 
predictions  agreed  well  with  the  observed  conditions.  ABL 
conditions  were  unstable  throughout  the  case.  Table  7  lists 
the  initialization  data  utilized  for  this  model  run,  and 
Table  VI  lists  the  2.5  hour  averaged  winds.  Pig.  4.7 
compares  the  two  input  winds. 


TABLE  7 

Case  Three  Initialization 


Atmospheric  Values 


Pressure: 
Pot.  T(C): 
Lift.  Cond 


1022.3 
17.7 
.  Lvl: 


mb 

508.0m 


Sea  Temp  (C)  : 
Spec.  Hum:  9. 
Aavection:  0. 


18.70 

48 


Mixed  Layer  Depth:  _ 

Below  Layer  Gradient:  -.00 


Ocean  Values 


11.5m 


5?mP. 


liie(E£T) 

1600 

1755 

2100 

2300 

0046 

0155 

0445 

0715 

0915 

1045 


Wj,nds 
Speed  { m/s) 

5.2 

6.5 

5.6 

4.5 

3.7 

5*2 

2.6 
2.  1 

4.8 
2.4 


Jump  (C) :  -.1 


Direction 

340.0 

338.0 

351.0 

352.0 

338.0 

325.0 

020.0 

296.0 

333.0 

271.0 


This  case  had  no  fog  related  problems  affecting  the 
dissipation  data  but  for  some  periods,  the  low  relative  wind 
speeds  degraded  dissipation  reliability.  Several  large 
dissipation  data  gaps  resulted.  For  the  bow  mast,  data  gaps 
exist  for  1600  to  0018,  0100  to  0815,  0845  to  1200,  and  from 
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Sp«ed  <M/S) 


Fig.  4.8  Case  III  Friction  Velocities 


1331  onward.  Only  seven  reliable  dissipation  values  were 
calculated  for  the  bow  mast.  Data  was  not  available  for  the 
main  mast  from  1600  to  0018,  0100  to  0815,  and  from  1410  to 
the  end  of  the  period.  This  resulted  in  four  thirty-minute 
dissipation  averages.  Fig.  4.8  compares  the  different  fric¬ 
tion  velocities. 

For  the  entire  period,  the  peak  and  averaged  wind  model 
u*'s  are  very  close  due  to  the  small  wind  differences 
between  the  two  (Fig.  4.7).  Both  the  main  mast  and  bow  mast 
dissipation  u*'s  agree  reasonably  with  u*'s  from  the  peaked 
wind  model  run.  This  is  particularly  true  in  the  second  half 
of  the  period.  Some  variation  is  expected  since  the  dissi¬ 
pation  u*  comes  from  instantaneous  measurements,  while  the 
model  uses  interpolated  values  between  the  input  winds. 
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Fig.  4.9  Case  III  Hixed  Layer  Depth 

The  model  also  predicted  the  MLD's  well  (Fig.  4.9) 
although  there  is  greater  deviation  from  the  observed 
(dashed  line)  than  in  Case  II.  As  before,  the  actual  HLD 
shows  some  internal  wave  or  advective  features  (particularly 
around  0000  and  0400).  Good  observed  MLD  data  was  available 
for  the  entire  period. 

The  initial  model  run  (solid  line)  agrees  well  in  the 
initial  stages  but  deviates  during  the  above  mentioned 
internal  wave  or  advective  events.  The  sharp  "dip”  at  about 
0930,  not  present  in  the  model  output,  is  an  artifact  intro¬ 
duced  by  use  of  a  spline  interpolation  in  the  plotting 
routine.  The  model  also  shallows  to  the  proper  depth  but  it 
occurs  about  one  hour  early. 
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Comparing  the  averaged  winds  MLD  with  the  peak  winds  MLD 
shows  that  both  are  very  similar.  The  plot  of  input  winds 
(Fig.  4.7)  shows  that  both  the  initial  and  averaged  winds 
were  very  close  throughout  the  period  except  between  0800 
and  1000,  so  both  model  runs  should  produce  similar  MID’s. 
After  1100,  the  averaged  wind  MID  (Fig.  4.9  dotted  line) 
shows  a  double  maxima  during  shallowing.  It  is  unknown  at 
this  time  why  the  model  predicted  such  maxima.  They  have 
been  slightly  exaggerated  by  use  of  the  spline  interpolation 
but  are  present  in  the  model  output.  As  with  the  initial 
wind  MLD,  an  anomolous  "dip"  is  observed  at  1030  also  due  to 
the  spline  interpolation.  Both  the  initial  and  averaged 
wind  MID's  are  within  operational  tolerances  for  Navy  use 
but  the  forecaster  should  be  well  versed  in  the  model  prop¬ 
erties  and  limitations. 
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7.  CONCiaSIONS  AND  RBCOHMENDATIOSS 


A.  CCJC1D5IOIS 

The  coupled  model  has  been  subject  to  various  past  anal¬ 
yses  and  performed  well.  This  study  has  examined  the 
performance  of  the  model  in  regimes  that  stretch  its  capa¬ 
bilities.  Several  important  conclusions  can  be  drawn: 

(1)  The  model  overpredicts  clouds  which  affects  short  and 
long  wave  radiative  effects  on  boundary  layer 
dynamics. 

(2)  Good  comparisons  of  bulk  and  dissipation  u*  for 
unstable  conditions  are  made. 

(3)  In  near  neutral  or  slightly  unstable  regimes,  u* 
predictions  vary  systematically  probably  due  to  swell 
and  wind  stress  interaction. 

(4)  Better  predictions  are  achieved  if  peak  winds  instead 
of  averaged  winds  are  used.  However,  a  forecaster 
would  only  have  average  wind  prediction  estimates 
available  at  the  time  of  operationally  initializing 
the  model. 

(5)  Given  reasonable  fluxes  cf  short  and  long  wave  radia¬ 
tion,  the  model  predicts  routinely  observed  varia¬ 
tions  in  ocean  mixed  layer  depth. 

(6)  It  is  a  potentially  powerful  tool  and  an  aid  in 
single-station  forecasting  to  a  Naval  Oceanographer 
who  is  aware  of  the  model  limitations. 

The  model  performed  very  well  within  its  limitations. 
However,  further  wock  needs  to  be  done  before  it  can  be 
implemented  in  an  operational  role. 
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B.  BECCMHENDATIONS 


In  order  to  facilitate  operational  usage,  several  areas 
should  be  explored.  Some  of  the  aforementioned  problems 
could  be  alleviated  by  attempting  the  following 
recoin  menda  tion  s : 

(1)  Continue  research  and  development  of  the  coupled 
model  with  direction  toward  operational 
i up le mentation. 

(2)  when  clouds  form,  the  model  assumes  them  to  extent 
from  horizon  tc  horizon.  This  is  unrealistic  and  a 
statistical  probability  of  percentage  cloud  cover 
should  be  added.  Partial  cloud  cover  would  also 
alter  short  and  long  wave  effects  on  boundary  layer 
dynamics. 

(3)  Incorporate  swell  and  wind  interaction  effects  into 
u*  calculations  based  on  bulk  formulae. 

(4)  Continue  to  analyze  HI1EEX  data,  particularly  for 
periods  where  there  are  more  complete  dissipation 
u*'s. 

(5)  Plan  for  model  coupling  to  other  predictive  models 
such  as  IBEPS  and  ICAPS. 

(6)  Involve  students  in  at-sea  data  gathering  experi¬ 
ments.  This  would  enhance  understanding  of  thesis 
research  and  minimize  "false  starts"  due  to 
attempting  cases  with  faulty  or  limited  data. 

The  model  has  tremendous  potential.  Performing  further 
analysis  and  tweaking  the  model  will  result  in  u  powerful 
tool  for  the  fleet;  one  that  can  be  used  in  times  of  crisis 
when  extensive  satellite  networks  may  not  be  available  for 
use  by  forecasters. 
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